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Abstract 

Thermal and electric conductivities are calculated for degenerate electrons scattered 
by phonons in a crystal made of atomic nuclei. The exact phonon spectrum and the 
Debye-Waller factor are taken into account. Monte Carlo calculations are performed for 
body-centered cubic (bcc) crystals made of C, O, Ne, Mg, Si, S, Ca, and Fe nuclei in the 
density range from 10^ to lO^i g cm at temperatures lower than the melting temperature 
but higher than the temperature at which the Umklapp processes begin to be "frozen out". 
A simplified method of calculation is proposed, which makes it possible to describe the 
results in terms of simple analytic expressions, to extend these expressions to any species 
of nucleus, and to consider face-centered cubic (fee) crystals. The kinetic coefficients are 
shown to depend tangibly on the lattice type. The results are applicable to studies of heat 
transfer and evolution of the magnetic field in the cores of white dwarfs and in the crusts 
of neutron stars. The thermal drift of the magnetic field in the crust of a neutron star is 
discussed. 
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1 INTRODUCTION 



Numerous studies have been devoted to calculations of the thermal and electric conduc- 
tivities of electrons scattered by phonons in dense crystalline matter of the cores of white 
dwarfs and the crusts of neutron stars (see, e.g., Yakovlev and Urpin 1980; Raikh and 
Yakovlev 1980; Itoh etal., 1984, 1993 and references therein). In particular, Raikh and 
Yakovlev (1982) calculated the thermal and electric conductivities by the Monte Carlo 
method taking into account the exact phonon spectrum in a Coulomb crystal of atomic 
nuclei, as well as the contribution of the Umklapp processes and the normal electron- 
phonon scattering processes. 

Later, Itoh etal. (1984, 1993) performed new calculations for the same conditions 
taking into consideration the Debye - Waller factor. This factor describes suppression 
of the electron-phonon scattering for large amplitudes of ion vibrations in crystals. In 
their calculations, the authors used an approximate treatment of the Umklapp processes. 
Moreover, they fitted their results by very cumbersome formulas. In addition, all the 
quoted calculations were made for bcc Coulomb crystals, whereas fee crystals may also 
exist in stellar matter (Section 6). 

The aim of this study is to repeat the exact ("from the first- principles") Monte Carlo 
calculations of Raikh and Yakovlev (1982) including the Debye - Waller factor, to study 
the dependence of the kinetic coefficients on the lattice type, and to present the results 
in the form suitable for practical use. Possible applications of the results are outlined, 
including the thermal drift of the magnetic field in the crust of a neutron star. 



We will restrict ourselves to the case of sufficiently dense matter, p > AZ g cm~^ (A 
and Z are the mass and charge numbers of nuclei, respectively), where the electron gas is 
degenerate and almost ideal, while the atoms are completely ionized by electron pressure. 
For simplicity, we consider only one atomic species. 

The state of degenerate electrons is described by the Fermi momentum pp = fik-p = 
^(Svr^ne)^/^ or by the relativistic parameter 



where rig is the electron number density, mg is the electron mass, pg is the density in 
units of 10^ g cm and is the number of baryons per electron. The electron gas is 
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nonrelativistic (x <C 1) for p <C 10® g cm ^ and becomes ultrarelativistic (x 3> 1) for 
p > 10® g cm"3. 

The state of ions (nuclei) is characterized by the parameter of nonideaUty V = 
Z^e^/{akBT), where a = [3/(47rrai)]^/^ is the mean separation between the nuclei, is 
the number density of the nuclei, and kB is the Boltzmann constant. At moderately high 
temperatures, the nuclei form a Coulomb crystal. A bcc crystal, which is studied in the 
main part of this paper (a fee crystal is analyzed in Section 6), is the most tightly bound 
crystal (sec, e.g., Brush etal. 1966). We will consider the temperature range Tu <T < T^. 
Here, Tu ~ TpZ^/^e^ / {hvp) is the temperature below which the Umklapp processes are 
frozen out (see below), fp is the Fermi velocity of the electrons, and and Tp are the 
melting and ion plasma temperatures, respectively, 

Tm = 4^«. 1.323 xl0^zV3(M'^'i!^K, (2) 

akBi-m KfJ-eJ i-m 

Tp = ^ « 7.832 X 10% /I^K. (3) 

kB V ^'"e 

The crystal melts at F = F^ ~ 172 (Nagara etal. 1987). Under conditions of study, 
Tm is significantly lower than the electron degeneracy temperature. Strictly speaking, 
the melting temperature of crystals made of the lightest nuclei (H, He) decreases at high 
densities due to strong zero-point vibrations of the nuclei (Mochkovitch and Hansen 1979). 
We will not consider these nuclei. The plasma temperature Tp is determined by the ion 
plasma frequency ujp = ^/AiTZ'^e'^ni/mi, where rui is the nuclear mass. Note that the Debye 
temperature of the crystal is Td ~ 0.45rp (Carr 1961). At T >Tu, the electron-phonon 
scattering can be treated in the approximation of almost free electrons. If Z » 1, the 
main contribution to the scattering comes from the Umklapp processes. The matter is a 
good conductor; there are many free electrons per one ion. The latter circumstance clearly 
distinguishes dense stellar matter from "terrestrial" metals (see, e.g., Yakovlev and Urpin 
1980). At T <^ Tjj, the character of the electron-phonon scattering changes (Raikh and 
Yakovlev 1982). The Umklapp processes are frozen out: as the temperature decreases, 
these processes occur only in those cases where the electron momenta before and after 
scattering lie near intersections of the electron Fermi surface with the boundaries of the 
Brillouin zones (where the electron energy spectrum contains an energy gap). We do not 
consider temperatures T -^Tu as they are usually not important for applications. 

The thermal k and electric a conductivities are conveniently expressed in terms of the 
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electron effective collision frequencies i^^ and 

K = « 4.04 X IO^^x^PTq ' ' ^ 



Sm^h'f^ \ i^K I cm s K 

, . . 149 X 10-.^;3 f 1. (4) 

Here, /3 = t^p/c = + x^, m* = me\/l + x^. 

To calculate z/^ and z^o-, we use the well-known variational method (Ziman 1962) with 
the simplest trial functions that describe deviation of the electron distribution from equi- 
librium. For astrophysical conditions, this method was used, for example, by Flowers and 
Itoh (1976) and all their successors. At T > Tp, this method yields an exact result, while 
at T ^ Tp, it gives an error of no more than a few percent in calculations of and no 
more than a few dozen percent in calculations of f^. 

We describe the electron states using the scheme of extended Brillouin zones in the 
approximation of free electrons. In this case, the Fermi surface and the dispersion relation 
take the same form as those for the free electron gas. We obtain 

= ^ ^i^<x,« ~ 0.955 X 10^5 s-\ (5) 

where we have introduced convenient dimensionless functions 

d5dS' 



7? 



t^S^ JsJs q^HqW 



g-2W(g) 



in which 



|/(g)P5:[ge.(fe)]2^^^^5<x,., (6) 



t = — = O.USTq , , , , 

9. = q , 9. = q + (7) 

Here Tq is the temperature in units of 10^ K, and s = 1, 2, 3 enumerates three branches 
of phonon vibrations. Integration in (P) is carried out over all possible positions of the 
momenta p and p' of an electron before and after scattering on the Fermi surface; S = 
Airkp'^ is the Fermi-surface area. The quantity hq = p — p' is the electron momentum 
transfer in a scattering event, and k and es{k) are, respectively, the wave vector and the 
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polarization unit vector of a phonon excited or absorbed by an electron: ±k = q — K, 
where K is such a reciprocal lattice vector that k lies in the first Brillouin zone. For k = q 
{K=0), here occur the so-called normal scattering processes, while for k ^ q {K ^ 0), one 
has the Umklapp processes. The quantity e{q) in (^) is the static longitudinal dielectric 
function of the electron gas. It describes the screening of the Coulomb potential of nuclei 
by electrons. The quantity f{q) is the nuclear form-factor that takes into account finite 
size of atomic nuclei. Finally, W{q) is the Debye-Waller factor (see, e.g., Davydov (1983)) 

-(..--(.)^^i:/(|J,J^-(|), „ 

where integration is carried out over the first Brillouin zone. Calculations show that under 
the conditions of interest, it is sufficient to set 

where Tq is the rms deviation of an ion at a crystal site. The angular brackets denote 
averaging over phonon frequencies and polarizations: 

{fsm = ^Y.I dfe/,(fe), (10) 

Here, Vb is the volume of the Brillouin zone. At t ^ 1, we have (2;~^cth {z/2)) = 2n_2i^) 
while at t < 1, (z^^cth (z/2)) = u-it. In this = {(io/ujp)'^) are the frequency 

moments of the phonon spectrum. 

According to (^) and (|9|), the Debye-Waller factor weakens the electron-phonon in- 
teraction and increases the thermal and electric conductivities. Clearly, it is important 
for strong ion vibrations in the lattice - near the melting point (when thermal vibrations 
are especially strong) and at high densities of the matter (when zero-point vibrations are 
strong) . 



3 METHOD OF CALCULATION 

Thermal and electric conductivities are calculated by 4D integration in (^) over all orien- 
tations of the electron momenta p and p' on the Fermi surface. The Monte Carlo method 
has been used for numerical integration. 

Initially we have calculated the phonon spectrum for a bcc Coulomb crystal immersed 
in a uniform compensating electron-charge background. The frequencies LOsik) and polar- 
ization unit vectors e<j(fc) of phonons of three modes {s = 1, 2, 3) are determined from the 
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system of equations composed of the elements of the dynamic matrix (Kohen and Keffer 
1955; Carr 1962). Extensive tables of the dynamic matrix elements have been calculated 
for values of k in different points of a primitive cell of the first Brillouin zone. To de- 
termine the phonon parameters for an arbitrary wave vector k, we have interpolated the 
tabulated data and then solved the system of equations. Near the center of the Brillouin 
zone, for k <^ q£) (where qd = (Ovr^nj)^/'^ is the radius of the sphere with a volume equal to 
that of the Brillouin zone), we have used the well-known asymptotic expressions (Kohen 
and Keffer 1955) for the dynamic matrix elements. In the approximation of a uniform 
electron background, for k <C gz), two phonon modes (s = 1 and 2) prove to be acoustic 
{us ~ ^pk/qo)-, while the third {s = 3) is optical {ujs ~ oJp)- Response of the electron 
background is important only for the frequency of the third mode for k <^ q^. Under 
the action of the response, the optical mode transforms into an acoustic one (Pollock 
and Hansen 1973). To allow for this effect, the frequency W3(fc) obtained for the uniform 



electron background has been multiplied by k/J k"^ + fexF^) where 



/ 4e2 

kTF = kFd^ (11) 
y TrnvF 

is the inverse radius of electron screening of a charge in plasma. 

In our calculations, we have used the static longitudinal dielectric function £]{q) of a 
degenerate relativistic electron gas obtained by Jancovici (1962). The simplest (see, e.g., 
Itoh etal. 1993) form- factor has been chosen, 

3 

/(?) = 7 [sm{qrc) - qrcCOs{qrc}] , (12) 

[qrcr 

corresponding to the uniform charge distribution in the proton core (of radius rc) of an 

atomic nucleus. In the density range of interest (see below), the radii of nuclei do not 

change under the ambient pressure, and we can use the standard formula rc = 1.15 x 
20-13^1/3 

Initially we have calculated the Debye - Waller factor (^ for the bcc lattice. The 
results closely match those of Itoh @et al.@ (1984, 1993). However, we have obtained a 
simpler fitting equation: 

, , q^fi ( , , 4.9801 + 13.00 X 61.099^2 \ ^ ^ 

2W{q) = 1.3995 exp -2.7256t -, ^ „qQ,2 • 13 

rriiUJp y 1 + 61.099t"^ J 

This expression gives an error less than 1% for any q. We have calculated the integrals (^) 
for eight elements: ^^c, i^q, 20^6, '^'^Mg, ^^Si, 223, ^OCa and ^^Fe. Unlike Itoh etal. (1984, 
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1993), we have not considered H and ^He. For hydrogen, the Umklapp processes are not 
allowed at all, and the kinetic coefficients require a separate treatment. The helium melting 
temperature (Mochkovitch and Hansen 1979) is likely to be lower than the temperature 
Tfj at which the Umklapp processes are frozen out, which is beyond the scope of our study 
(Section 2). We have calculated and Fo-F@k for the electron relativistic parameters x 
= 0.1, 0.3, 1, 3, 10, and 30 and for the dimensionless temperatures t = 0.02, 0.04, 0.08, 
0.2, 0.4, 0.8, 2, 4, 8, and 20. The selected values of x correspond to the density range from 
approximately 2 x 10^ g cm" 3. Unhke Itoh etal. (1984, 1993), we have not considered 
lower densities p, since the approximation of complete ionization does not hold at these 
p (see, e.g., Yakovlev and Urpin 1980). For carbon, we have excluded the value x = 30, 
because at such a high density carbon is instantaneously burned in pycnonuclear reactions 
(Yakovlev 1994). For all elements, we have also excluded those t at which T > 2Tm (see 
(^)). Calculations of and F^- for each element at fixed x and t have required up to 
500 000 configurations in "Monte Carlo chains" (choices of orientations of the momenta p 
and p' on the Fermi surface) to achieve an acceptable accuracy of 5 - lOIn all the cases of 
interest, an inclusion of the nuclear form- factor (allowance for the non-point distribution 
of the nuclear charge) has virtually no effect on the final result. The finite sizes of nuclei 
become important only at p > 10^^ g cm~3. 



4 NUMERICAL RESULTS AND THEIR ANALYTIC ANAL- 
YSIS 

For practical applications, it is useful to fit the functions F„ and F^ by simple analytic ex- 
pressions similar to those proposed by Yakovlev and Urpin (1980) and Raikh and Yakovlev 
(1982). Under the conditions of interest, the main contribution to the electron scattering 
comes from the Umklapp processes. The Fermi surface is crossed by many boundaries of 
the Brillouin zones, and we may set in (^) 

El««'(*)l'(ir^ " <i=«"t"o,.(*), (14) 

where n = or 2, and the angular brackets imply averaging (|10|). 

The functions Gn{t) were calculated by Yakovlev and Urpin (1980), who also proposed 
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simple fit formulas: 



Go{t) 



ao 



The coefficients in (^) are chosen in such a way to reproduce the asymptotic expressions 
for t <C 1 and t ^ 1; equations (|l6|) describe G„(t) with an error less than lOfrequency 
moments of the phonon spectrum (Section 2). For a bcc lattice, we have uL^a'''' = 13.00 
(see, e.g., Mochkovitch and Hansen 1979), af''''^ = 0.0174 and a^^"""^ = 0.0118. 
Using the approximation (|T4[), from (^) we obtain 

Following Yakovlev and Urpin (1980), we have taken into account that the simplified 
expression (14) is valid only for the Umklapp processes. Therefore, the region of integration 
over the electron momentum transfer g in a scattering event is limited by a minimum 
momentum qmin- 

We set Qmin = = (evr^nj)^/^, where go is the radius of a sphere with 
the volume Vb- For q < Qb, the normal scattering processes dominate. 

ID integrals (|l^) are readily calculated numerically. However, they can also be approx- 
imately expressed in an analytic form. Indeed, finite sizes of nuclei are unimportant, under 
the conditions of interest (Section 3), and we may take /(g) = 1. The dielectric function 
e(g) of an electron gas differs slightly from unity for q > /ctf (pT|). Only for q <^ k^p does 
the quantity |e((/)P sharply increase: the electron screening becomes important, giving 
rise to an effective cutoff of the integrals for low q < k^p. Simple estimates show that 
/ctf ^ in all the cases of study except for those of the lowest density {p ~ 10^ g cm^^). 
At p < ZA g cm^^, the electron screening would become extremely important, but our 
approach becomes invalid. In particular, the condition of complete ionization is violated 
(Section 2). In calculating the integrals (^) approximately, the electron screening may be 
taken into account by shifting the lower integration limit qmin- We set q^-^ ~ + ^tf^ 
and e = 1 for q > qmin- The integrals can then be taken: 

F^ = Go{t)[2Ro{s^)-l3^Ri{s^)], (18) 

F, = Go{t)[2Ro{s,)-(3^Ri{s,)] 

+ ^G2it)[f3^Riis^)-3(3''Ro{s^)-2Rois^) + 3R2isi^)]- (19) 
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Here, we have introduced the functions 

Rois) = -(e-°^-e 
a 

2 



-a\ 



= [e-"'^(l + as)-e-"(l + a)], 



a 



R2{s) = E{as)-E{a), (20) 

in which E(x) is the integral exponent, a = Ar^kp'^ is defined by in the Debye-Waher 
factor(^ , and s = q'^j„/(2A;F)^. It is convenient to take somewhat different lower limits 
qmin in the functions and F^- (see below). In our case, qmin kp and s <^ 1. Note 
that E((7) can be calculated conveniently using the simple fit formula 



OO 

y 



exp 



g3 + 0.1397 



1\ 0.5772 

In l + 



1 + 2.2757g2 



(21) 



With a mean error of about 1%, the Debye-Waller parameter a can be fitted by 



a = ao [lu-ie-^-^^^^ + tu-2] , ao = ~ 1-683J-|^, (22) 



2 / kp,Tpmi V 

where U-\ is another frequency moment of the phonon spectrum. For a bcc lattice, 
^(6cc) ^ 2.800 (Mochkovitch and Hansen 1979). 

If a <C 1, the Debye-Waller factor is unimportant. Since s <C 1, in this case we obtain 

i?o(s)«i?i(s)^l, i?2~lnQ). (23) 

These formulas correspond to the approximation used by Yakovlev and Urpin (1980), as 
well as by Raikh and Yakovlev (1982). The inclusion of the Debye-Waller factor weakens 
the electron scattering, and the functions i?Oi -^i and become smaller than (23). The 
functions ( pOD completely determine the dependence of electric and thermal conductivities 
on the Debye-Waller factor. 



Let us compare the approximate analytic expressions (18) - (^2|) for Ffj and with 
the results of our numerical calculations (Section 3). The values of appear to be fitted 
by equation (|18|) with an error not larger than the error of the computations if we set in 
the functions Ro{Sf^) and Ri{sa) 

Sa = SD + STF, (24) 
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If one uses the same parameter @s for thermal conductivity 

Sk = Sift = Sa = SD + STF, (27) 

formulas ( p!9|) - fit the calculated values of with a mean error of 10% for all eight 
chemical elements and all x and t. The approximation ( ]27| ) is least accurate for light 
elements, primarily for carbon (for it, the mean error is about 15%, and the maximum 
error of ~30% is obtained at x = 10 and t = 0.2). The accuracy of the Ff^ fit can be 
improved, on average, by a factor of about 1.5, if we take 

Sk = sdC + stfL 
siK = sdCi + stfCi, 

834 

C = 0.101 + 0.0305Z, ^ = 3.84 + — 5-, 

2.66 257.4 4874 , , 

Ci = ^0122 -1' ^1 = 1-456-^ + ^. (28) 



For applications, we recommend to use (p^) and (27) because their accuracy is quite 
acceptable. It is clear from the derivation of equations ( |2^ and (^) that they can also 
be applied to other chemical elements. Equations ( [2^ ) are more accurate, but they hold 
only for 6 < Z < 26. 

The functions and were previously calculated with allowance for the Debye- 
Waller factor by Itoh etal. (1984, 1993). These authors started from approximate equa- 
tions ([TtI) , whereas we have performed our numerical calculations using the general but 
much more complicated equations (w. Our calculations confirm the validity of simplified 



equations (17). Hence, our numerical results agree, within the calculational errors, with 
those of Itoh etal. (1984, 1993). We have thus proved the validity of the approximate 
calculations of electric and thermal conductivities based on equations ( |l^ ) and (|T7|). Fur- 
thermore, we have carried out an analytic analysis of equations (|l^). This has enabled us 
to obtain the simple fit equations (|T8|) - (^) valid to all chemical elements considered. 
As follows from their derivation, these equations are also valid for other elements (see 
above). Itoh etal. (1984, 1993) also fitted their results by analytic expressions, but their 
expressions are very cumbersome. They include 56 adjustable parameters for each element 
and do not allow one to extend the results directly to other elements. 
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The quantities Fcr and F^, calculated with and without allowance for the Debye- 
Waller factor, are shown in Figs. 1 and 2. Including the Debye- Waller factor increases 
the electric and thermal conductivities of matter. The effect is especially pronounced in 
the cases where the amplitude of ion vibrations in a lattice cite is large, that is, near 
the melting point (where thermal vibrations are strong) and at high densities (where zero 
vibrations are important; Fig. 2). At T ~ Tm the Debye- Waller factor increases the 
kinetic coefficients by a factor of 2 - 4. 

At t > 1, the electron-phonon scattering becomes quasi-elastic (the energy transfer in 
a scattering event is much smaller than kBT). Using ([ITD, ( [20| ) and ([2^), one can easily 
show that in this case a single effective collision frequency can be introduced for both 
electric and thermal conductivities: 

kBT 

2u 



a(l-e-") (!-(! + a) e-°) . (29) 



We have taken into account that sa <C 1 in virtually all the cases studied. Equation ( p9| ) is 
a generalization of the similar expression obtained by Yakovlev and Urpin (1980) without 
allowance for the Debye- Waller factor. 

For t ^ 1 and q ^ 1, equations (]2^ ) and (p9|) yield the simple expression 

In this limiting case, the collision frequency of electrons is temperature independent and 
close to a similar frequency (see, e.g., Yakovlev and Urpin 1980) in the Coulomb liquid 
(T > Tm). Formally, the collision frequency (|30| ) would be equal to the frequency in the 
liquid if the Coulomb logarithm were equal to A = 0.5. 



^0 = ^7^. (30) 



5 THERMAL DRIFT OF THE MAGNETIC FIELD 

The effective collision frequency (|29| ) defines, in particular, the thermal drift of the mag- 
netic field B in the neutron star crust (Urpin and Yakovlev 1980). Recall that the ther- 
mal drift is caused by the heat flux emerging from the stellar interior. For simplicity, 
we consider a weakly magnetized degenerate electron gas in which cobTo *C 1. Here, 
lub = eB/{m^c) is the gyrofrequency of electrons with energies equal to the Fermi energy, 
and To = l/i^o is the effective relaxation time. Under these conditions, the heat flux is 
almost unaffected by the magnetic fleld, and the thermal drift is determined by the Hall 
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component of the specific thermal electromotive force. The thermal drift velocity u is 
parallel to the heat flux Q and for t > 1 is equal to 



u 



V 



2r]Q 



72r] 



UePFVF 

2 



Vl06 K 



m 

yr' 



dlnp 



(31) 
(32) 



Here, we set Q = uTg , where is the effective temperature of the stellar surface. The 
numerical coefficient ~ 1 is determined by the logarithmic derivative with respect to 
the electron momentum. Once the derivative in (32) is taken, p = pf should be set. 
Basically, rj can have different signs, depending on the mechanism of electron scattering. 
In the Coulomb liquid (T > T^), according to Urpin and Yakovlev (1980), rj > 0, and the 
thermal drift is directed outward. In crystalline matter, from (|32|) at t > 1 we obtain 



V 



1 



a 
3 

^"2 = -3 



(1 



n 



2«/3V2 + 3(/3' - 

[l-(l + «)e-]. 



6ro 
2 



(2 + 2Q + a^)e" 



(33) 



(34) 



Plots of Tj against density in the neutron star envelope with the temperature T = 8 x 10^ 
K are shown in Fig. 3. This temperature corresponds to Tg ~ 10^ K. The coefficient r] 
calculated for the crystal without allowance for the Debye- Waller factor is indicated by 
long dashes. It can be obtained from and (|3|) in the hmit a ^ 0: r/ = (2 - 3(3^) /{A - 
2/32 ) (Urpin and Yakovlev 1980). In this approximation, r] is positive in a low-density 
electron gas (for Vp/c^ < 2/3) but becomes negative at higher densities. However, Fig. 
3 shows that at T ~ with allowance for the Debye-Waller factor, ij remains positive 
and close to a similar coefficient in the liquid phase. With increasing density and at a 
fixed temperature, t and r] decrease. Formally, at a sufficiently high density, equation (32) 
gives negative ry, but this occurs only at t < 1, when the formula becomes invalid. In the 
validity range of (^) and for typical values of the parameters, r] remains positive, and the 
thermal drift is directed outward. The calculation of at t < 1 is very complicated and it 
is beyond the scope of this study. 
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6 FACE-CENTERED COULOMB CRYSTALS 



The preceding astrophysical calculations of electric and thermal conductivities have been 
aU performed for body-centered Coulomb crystals that are most tightly bound. In addition 
to bcc crystals, fee and simple cubic Coulomb crystals can in principle exist. The binding 
energy of an ion in a Coulomb crystal of any type is U = —C^Z'^e^ja (Brush etal. 1966), 
where a is the ion sphere radius (Section 2), and the coefficients C are equal Cf)cc= 0.895929, 
C/cc=0. 895874 and Csc=0. 880059. The binding energies for bcc and fee crystals are very 
close to each other. The difference in the ion energies (in units of thermal energy) is 
(C^ftcc — Ufcc)/ikBT) ~ — 0.000055r, where F is the nonideality parameter (Section 2). 
Near the melting point, T T^, we have T « 172 and {Ubcc — U fee) / {k^T) ~ —0.01. 
Consequently, in addition to bcc crystals, fee crystals can be be formed in many cases. 

In Section 4, we have justified the validity of a simplified calculation of thermal and 
electric conductivities based on equations (p^) and (jl^). Let us apply this approach to 
fee crystals. Using ( p!5| ) and (p!7|), we can calculate the functions Fg- and F^. In this 
case, the type of lattice affects the functions G'o(t) and G2{t), as well as the Debye- Waller 
parameter a. Using the results of Kohen and Keffer (1955), we have created a computer 
code for calculating the phonon spectrum. The spectrum resembles the phonon spectrum 
of a bcc crystal. In particular, there are two acoustic and one optical modes near the 
center of the Brillouin zone, but the frequencies of acoustic phonons are somewhat lower. 

( fee) ( fee) 

As a result, the frequency moments u_i = 4.03 and u_2 = 28.8 are larger than those 

f fee) 

for a bcc crystal. We have also calculated some other frequency moments: =0.462, 
^(^/cc)_Q 2gy_ The integration over the Brillouin zone in the expressions for Go{t), G2{t) 
and a has been carried out by the simplified method proposed by Mochkovitch and Hansen 
(1979). The results can be fitted by the same equations (^) and (|2^ ) as for the bcc lattice 

( fee) ( fee) 

but with other values of constants: Cg =0.00505, =0.00461. The fit accuracy turns 
out to be of the same order of magnitude. All equations (|l8|) - ( p8| ) for F^ and Fo- remain 
valid. 

The plots of Ff^ and Fu versus t for the bcc and fee lattices made of iron nuclei are 
given in Fig. 4. The curves are similar but the peaks at t ~ 1 are much more pronounced 
for the fee lattice, which is a result of a softer spectrum of acoustic phonons. At higher t, 
the curves approach each other since they are described by the common asymptote (|30|) . 

A simple cubic lattice is appreciably more weakly bound than a bcc or fee lattice. 
In addition, the phonon frequencies in some crystallographic directions of a simple cubic 
lattice turn out to be complex, i.e., such a lattice is unstable. This is a common and well- 
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known (Born 1940) property of cubic lattices composed of particles interacting through 
central forces. 

7 CONCLUSION 

We have calculated the electric and thermal conductivities of degenerate relativistic elec- 
trons scattered by phonons in Coulomb bcc and fee crystals made of atomic nuclei taking 
into account the Debye - Waller factor. The results of these calculations are valid at den- 
sities 10^ g cm~^ ^ P ^ 10^^ g cm~'^. For bcc crystals, the calculations are performed by 
the Monte Carlo method using the exact formulas. Our results confirm the validity of the 
approximate approach used by Itoh etal. (1984, 1993). An analytic analysis of the results 
has enabled us to describe them by the simple fit expressions ( p^ - (^) and extend the 
calculations to the fee lattice. It is shown that thermal and electric conductivities are 
sensitive to the lattice type. 

One might expect our results to be also applicable at sufficiently low densities, 1 g 
cm^'^ < p < 10"^ g cm^'^, where the pressure ionization of atoms is incomplete. In this 
density range, the charge number Z of nuclei should be replaced in all the above expressions 
with an effective charge number .^e// ^ Z determined by electron screening of the nuclear 
charge. However, no reliable calculations of Z^ff have been made as yet. 

Our results are valid for degenerate cores of white dwarfs and to outer degenerate 
envelopes of neutron stars. First of all, these results can be useful in studying the thermal 
evolution of these objects (cooling, nuclear burning of accreted matter) and the evolu- 
tion of their magnetic field (ohmic dissipation, generation under thermomagnetic effects). 
In particular, we have shown in Section 5 that an inclusion of the Debye-Waller factor 
strongly affects the thermal drift of the magnetic field in the crust of a neutron star. 

The results can easily be generalized to higher densities, 10^^ g cm^^ ^ P ^ 10^"^ g 
cm~^. To do this, the functions F^^ and should be calculated with allowance for finite 
sizes of atomic nuclei for realistic models of superdense matter of subnuclear density. We 
plan to consider this in a separate work. 
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Figure 1: Fo- (curves 1) and (curves 2) versus t = T/Tp for a crystal composed of ^^Fe 
at the density p = 2.1 x 10^ g cm~^ (x = 1) with (sohd hues) and without (dashed hues) 
allowance for the Debye-Waller factor. 
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Figure 2: Same as in Fig. 1 but for ^^C nuclei at p = 5.26 x 10^ g cm~^ {x = 3). 
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Figure 3: The coefficient rj that determines the thermal drift velocity of the magnetic 
field (see (^l])) as a function of density at T = 8 x 10^ K for matter composed of ^^Fe 
nuclei. Curve 1 corresponds to the Coulomb liquid (Urpin and Yakovlev 1980), curve 2 
to a crystal without allowance for the Debye-Waller factor, and curve 3 to a crystal with 
allowance for this factor. The vertical line shows the density above which solidification 
occurs. 
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Figure 4: Same as in Fig. 1, for bcc (solid lines) and fee (dashed lines) crystals with 
allowance for the Debye-Waller factor. 
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